Multiple fluxoid transitions in mesoscopic superconducting rings 
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The authors report magnetic measurements of fluxoid transitions in mesoscopic, superconducting 
aluminum rings. The transitions are induced by applying a flux to the ring so that the induced 
supercurrent approaches the critical current. In a temperature range near T c , only a single fluxoid 
enters or leaves at a time, leading to a final state above the ground state. Upon lowering the tem- 
perature, several fluxoids enter or leave at once, and the final state approaches the ground state, 
which can be reached below approximately 0.5 T c . A model based on the widely used time depen- 
dent Ginzburg-Landau theory for gapless superconductors can only explain the data if unphysical 
parameters are used. Heating and quasiparticle diffusion may be important for a quantitative un- 
derstanding of this experiment, which could provide a model system for studying the nonlinear 
dynamics of superconductors far from equilibrium. 

PACS numbers: 74.40.+k, 74.78.Na 



The dynamics of phase slips (PS) in quasi-lD su- 
perconducting wires has long been studied intensively. 
Most observed phenomena have been explained with 
phenomenological models [3, 0| and the time-dependent 
Ginzburg-Landau theory (TDGL) for gapped supercon- 
ductors [3|, often with good quantitative agreement 
However, in those wires with diameter d of typically a 
few /Ltm, the quasi-lD limit £(T) > d and A(T) 3> d was 
only accessible very close to the critical temperature T c , 
where the coherence length £ and penetration depth A 
diverge. Due to various approximations, the microscopic 
validity of most of the theoretical models employed is 
also limited to the vicinity of T c , Modern e-beam and 
nanotcmplating techniques yield superconducting wires 
which remain quasi-lD down to T = 0, and can be thin 
enough (d ~ 10 nm) to study phase transitions due to 
quantum PS's. Such quantum phenomena, which have 
recently attracted tremendous interest [f|, are possible 
when some action S ~ AEAt is comparable to h. Thus, 
the dynamical laws determining the time scale At are 
very important for understanding quantum effects. Fur- 
thermore, a better understanding of superconducting dy- 
namics is interesting both for its own sake and for device 
applications, such as microbridgc SQUIDs 0] or nanowire 
single photon detectors Q- 

We have probed the dynamical evolution of the or- 
der parameter in mesoscopic Al rings at temperatures 
T = 0.15 . . .0.9 T c . Upon increasing the magnetic flux 
$ a threading a ring, the circulating supercurrent I in- 
creases quasistatically until it approaches the critical cur- 
rent, and the ring switches rapidly to a more stable state. 
We find that near T c , the phase winding number n of the 
order parameter ip(x) = |?/>|e m:r changes by An = 1 at 
each transition, which can leave the ring in a metastable 
state well above the ground state. At lower T, An in- 
creases, and the final state approaches the ground state, 
while static states persist over a larger range of applied 



flux. Below T ss 0.5 T c , the ring decays into the ground 
state over an increasing range of flux bias points. Such 
multiple fluxoid transitions were predicted based on the 
gapless TDGL equations [§] and observed at a single T 
Q. Refs. [lj| and 11 1 contain further related theoreti- 
cal work on state selection and the effect of fluctuations 
on the dynamics of mesoscopic superconducting rings, al- 
beit in a different parameter range. The TDGL model 
in [1, 0] can only reproduce our and previous ring data 
by choosing a numerical constant which is inconsistent 
with microscopic theory and transport measurements in 
the resistive state. The gapped TDGL theory may over- 
come this inconsistency. However in Al, both versions of 
TDGL are justified only extremely close to T c , and we es- 
timate that the neglected heating and quasiparticle (QP) 
diffusion may be important in our experiment. Lacking a 
tractable theory for superconducting dynamics far from 
T c , it is common to use TDGL, sometimes well below T c 
and in the gapless version. Thus, it is of considerable 
interest to understand to what extent those theories are 
adequate outside their theoretical range of validity, and 
to develop better theories. Experiments such as ours can 
be used for validation. 

Unlike transport measurements on superconducting 
wires, where voltage reflects the mean phase wind- 
ing rate, our experiment probes when phase unwinding 
ceases, once started. Thus, we explore fast dynamics 
without a high measurement bandwidth. Furthermore, 
DC transport measurements are prone to thermal run- 
away effects at low T, where the critical current is large 
and cooling is inefficient. In our case, heating is also sig- 
nificant, but localized in space and time due to the finite 
number of PS's. The absence of leads eliminates the need 
to consider proximity effects and boundary conditions at 
the contacts. 

The results reported here are from 13 rings with radii 
R = 1 and 2 /nm and line width w = 70, 115, 140 and 



2 



180 nm, measured in a single cool down. We measured 
7 as a function of $ a and T by positioning a scanning 
SQUID microscope [12J over each ring individually. The 
samples were made by e-beam lithography using PMMA 
resist and liftoff. The 60 nm thick Al film was deposited 
using e-beam evaporation at a pressure of 10 -7 mBar and 
a rate of 35 A/s. Using £ = 1-6 /im for pure bulk Al [l3} 
and £(0) = 140 to 180 nm depending on w, obtained as 
described below [see also Fig. QJg)], we infer a mean free 
path of l e = 1.4£(0) 2 /£ = 17 to 28 nm. The measured 
T c is 1.25 K. 

Our SQUID sensor is designed as a susceptometer as 
discussed in 14. We positioned one of its pickup 

loops over a ring and recorded its response while sinu- 
soidally varying $ a , generated by an on-chip field coil, at 
0.6 Hz. To isolate the flux generated by I, we subtracted 
a background measured by retracting the SQUID from 
the sample. For each ring, we took 40 such <& a - 1 curves 
at each T, in steps of 10 mK — a total of more than 
40,000 field sweeps containing 550,000 transitions. 

We also collected data at a smaller $ a -amplitude, so 
that fluxoid transitions only occurred near T c , to extract 
\- 2 {T) and £(T) [Fig Eft g)], as described in Ref 0. 
Using fifth order rather than cubic polynomials to fit the 
low T <& a -I curves accounts for corrections to ID GL 
at low T0, and for imperfections in the rings, whose 
effect is biggest at large I. We also corrected for partial 
screening of the applied flux. 

Fig. Q] shows a few individual field sweeps. Each 
branch of a 3> a -I curve represents the response I„ (ip = 
^a/^o) °f a state with fluxoid number n, with a linear 
response around <p = n, and a curvature caused by pair 
breaking at larger ip—n. <&o = h/2e is the superconduct- 
ing flux quantum. Since w <C R and H a <C H C 2, I n and 
the transition points only depend on (p — n to high accu- 
racy. In the GL regime, I n {<p) = ~(wd&o/2irRfio\ 2 )((p — 
n)(l— (t; 2 /R 2 )(<p — n) 2 ). At lower T, deviations from the 
cubic form become noticeable. 

We have extracted the position tp t of each transition 
and n for each branch, and computed </>t(T) = (\ip t — 
by averaging over all transitions in each set of field 
sweeps, rii and rif are the n before and after each tran- 
sition, as illustrated in Fig. Eld). The < 50 m$o rms- 
scatter of \ipt~ n i\ around <pt, which is partly due to sensor 
vibrations, is negligible for our purpose. The observed 
frequencies of An = |rij — n/| values yield a probability 
distribution as a function of T for each ring, as shown in 
Fig- EI e )- Since we find at most two consecutive inte- 
ger values of An at any given T, the distribution is fully 
characterized by its mean (An)(T) [Fig. [2]. 

A GL stability analysis neglecting screening effects 
and fluctuations predicts <p t — \J R 2 j £q L + 1/2/ V3 ~ 
R/V3£(T) In practice, 4>t as shown in Fig. Ovaries 
significantly between nominally identical rings, and is up 
to 30 % smaller than the measured value of R/y/3^(T), 
with the largest deviation for the smallest w. The varia- 




p(An) 



cp = <j> /o 

T a 
100 



0.5 













50 


CD 
II 


ii 


II 


co 
ii 


C\l| 

ii 1 ii 







< 


< 

I 


1 


<i < 


1 

0.5 







0.2 


0.6 




1 




T(K) 



T(K) 



FIG. 1: (Color online) (a)-(d) Typical single-sweep $ a -I 
curves for a R = 2^tm, w — 180 nm ring. The ring current I 
is proportional to the SQUID response $ squid, with mutual 
inductance M coup ~ 1.2 $o/mA for this ring size. At 1.1 
K (a), An = 1 always. Around 0.96 K (b) and 0.95 K (c), 
An = 2 transitions become increasingly likely. At 0.28 K (d), 
An = 6, and the ring decays close to the ground state at 
each transition. The arrows indicate the field sweep direction 
and illustrate An, ru, n/ and 4>t for the transition at tpt = 
—9.7. (e) Probabilities for a change An in the fluxoid number. 
(f),(g) A -2 and £ as a function of T. 



tions suggests that imperfections contribute significantly 
to this discrepancy, but a deviation from GL theory far 
below T c and thermal activation or tunneling may also 
play a role. The temperatures at which An changes also 
show significant scatter. However, they correlate strongly 
with the value of <pt at that point for similar T. Conse- 
quently, the curves of <j>t(T) — (An)(T) = ±(ip t — rif) 
from all rings [Fig. [2] approximately collapse into a 
band whose width of « 1 is largely due to the discrete- 
ness of (An). Below 0.5 T c , <pt — (An) < 1/2 can occur, 
which implies that the final state is the ground state with 
\(pt — nf \ < 1/2. The lower T, the larger is the range of <pt 
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We now turn to discussing our results in terms of the 
TDGL equations for gapped superconductors: 



FIG. 2: (Color online) (a) Mean transition point <j>t (lines 
only) and (An) (lines with dots, each dot representing one 
dataset of 40 field sweeps) . The inset shows the same data for 
just the ring from Fig. [T]and a R = 1 (im ring, (b) (f>t - (An) 
as a function of temperature for all 13 rings. The three shaded 
regions indicate the difference |n/ — n groun d\ between the n of 
the final and the ground state. Since <j>t — (An) = ±(v?t — rif), 



it encodes both |n/- 



id\, and the dependence of the final 



state on how far from an integer flux bias point a transition 
occurs. Comparison with the inset shows that <f>t — (An) 
depends directly on the temperature rather than R/£(T). 



over which the ground state is reached. The absence of a 
clear w dependence of 4>t — (An) vs. T indicates that the 
effect of the self inductance and the thermal activation 
barrier, both of which are proportional to wd, is small. 

The results show no significant dependence on sweep 
frequency from 0.6 to 9.6 Hz. One concern is heating 
from the ~5 GHz, ^5 m$o Josephson oscillation applied 
to the rings by the SQUID. Varying the pickup loop-ring 
inductance M coup by a factor 5 by retracting the SQUID 
about 2 /jm had no significant effect. However, when the 
sensor chip touched the sample substrate, 4>t{T) flattened 
slightly below 0.35 K, changing by 0.1 at 0.2 K. Since this 
effect is negligible, the data was acquired with the SQUID 
touching the sample to avoid variations of M coup due to 
scanner drift and vibrations. 
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v/TTrW \ dt ' " ' 2 dt 

= (V-iA) 2 V + (l- M 2 )V (la) 
V 2 U = V-Re(^*(-iV + A)V>), (lb) 

where 7 = 2teAo(T)/S. A (T) is the equilibrium gap, V 
the electrostatic potential, and te the inelastic electron- 
phonon scattering time. Without magnetic impurities, 
theory predicts u = 7r 4 /14<(3) = 5.79 Q. Eqs. Q} are 
written in dimensionless variables, as defined in Ref. [l7l ]. 
We only note that the unit of time is tql = 2ksTh/ A 2 ,, 
and the dimensionless ring circumference is 2ttR/£(T). 
The latter might lead to an indirect T-dependence. How- 
ever, Fig. EJ^b) shows only a weak i?-dcpcndence of 
the range of possible <f>t{T) — (An)(T), whereas plotting 
(f) t —(An) against i?/£(T) (inset) gives two distinct bands. 
Thus, the T-dependence of 7 is more important than that 
of R/£(T), to the extent that Eqs. (p} are applicable. 

The mechanism underlying multiple PS's in meso- 
scopic rings was analyzed in Refs. 0, 0| , based on numer- 
ical solutions of Eqs. Q]with 7 = 0. Ref. Q identifies two 
timcscales: the duration of a single PS, r^, and the relax- 
ation time T|^,| of I ip\ ■ If Tct> ^ , the phase may unwind 
by several 2ir before \ip\ can recover after the first PS. 
Further analysis for wires with finite 7 UJ shows that 
is related to the voltage drop AU = IRps across the PS 
center via the Josephson relation: = $o/AV, where 
Rps is the resistance of the region around the PS over 
which the electric field decays and the QP current is con- 
verted to supercurrent. The extent of this region is given 
by the charge imbalance length Aq» oc y/j/u^(T). One 
obtains oc \J u/jtgl and tui oc ^utql- The above 
expressions are valid for 7 » 1. For 7 < 1, one should 
replace 7 with 1. 

If treating u as a phcnomcnological adjustable param- 
eter and setting 7 = 0, one thus finds that t^/tui de- 
creases with increasing u, so that a larger u leads to 
larger An Q. Reasonable agreement with aprevious 
experiment at 0.4 K was obtained for u = 48 [9j. How- 
ever, u 3> 6 is inconsistent with microscopic theory and 
experiments in the resistive state of quasi-lD supercon- 
ducting wires 0, Qi| • Those show that except very close 
to T c , £ < Aq» oc \T — Tel" 1 / 4 . Some experiments even 
found a constant Aq* . When using Eqs. {1} with 7 = 0, 
this implies u < 1 and that u decreases with decreasing 
T |. For example, Eqs. JTJ with u = 0.01 and 7 = 
have been used to model the resistive state (rij . 

This inconsistency in the effective u arises because 
7 is negligible only in the gapless limit, where phonon 
(or magnetic impurity) induced pair breaking causes a 
fast reaction of the QP population to order parame- 
ter variations. For pure Al, this limit only applies for 
T > T c 
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K. An increase of 7 oc y/\T — T c \ with 
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decreasing T on the other hand both reduces t^/tui 
and increases Aq*/^(T), in qualitative agreement with 
both ring and resistive-state experiments. However, Eqs. 
((T|) neglect heating and QP diffusion, which is valid for 
T 4>, T \ip\ "C te- Thus, even accounting for the slow charge 
imbalance relaxation by allowing 7^0 only extends the 
theoretical validity to T c - T < 0.1 mK for Al Q. 

The observation that ipt — rif is nearly independent of 
R implies that the current after a transition is approx- 
imately inversely proportional to R. For wires carrying 
a quasi-dc bias current I on the other hand, phase slip- 
ping stops once I drops below some limit at which 
7> = <S> a /R PS Ici ~ tim [13]. Within the TDGL pic- 



ture discussed above, this finite length effect could be 
explained if the normal-like length determining Rps is 
set by circumference rather than Aq* , which is formally 
plausible since 2Aq» > 2irR at T = 0. However, QP's 
can only diffuse a distance wDt$ -C tR during a single 
PS, so that it seems questionable if the whole ring can 
contribute to Rps- Alternatively, the local reduction of 
the critical current at the PS center could increase as 
more PS's occur, so that larger rings holding more flux- 
oids have a smaller final current. Well below T c , where 
Ao ^> fcfiT, the electronic heat capacity c(T) is very 
small, while the kinetic energy density is large. There- 
fore, dissipating the latter in a PS will lead to significant 
heating, i.e. excitation of the energy mode of the QP 
population, which increases for larger R. Thcrmalization 
with the lattice occurs on the scale of te > 50 ns [l9| . 
whereas assuming that Rp$ is at least that of a ring sec- 
tion of length £ gives < 10 ps. Thus, the electrons 
are approximately a closed system. If the kinetic en- 
ergy was converted to homogeneous heating, the resulting 
electronic temperature would be given, very roughly, by 
max(T, 0.8 K), because c(T) is small for T < 0.5T C . At 
this temperature, 4>t is still rather large, so that uniform 
heating would lead to smaller An than observed. How- 
ever, given the short diffusion length \JDt^ ~ £ <C 2irR, 
hot QP's will remain localized so that diffusive cooling 
may affect tui. 

A general description of time dependent superconduc- 
tivity, similar to the Usadel equations for the static case, 
was derived in Ref. [3]. Unfortunately, the resulting 
system of equations for A, the Green's functions and 
electron distribution is rather complex. Approximate so- 
lutions for specific cases were studied in Ref. [2(| [2l| . 
Similar computations might allow a direct comparison 
between our data and microscopic theory. A much sim- 
pler first step could be a similar study as in Q using Eqs. 
(JTJ) , but with 7 rather than u as free parameter. Heating 
and QP diffusion could be accounted for through a local 
effective temperature governed by the heat equation. 

In conclusion, we have shown experimentally that the 
final state after a field induced fluxoid transition in meso- 
scopic rings approaches the ground state for T — > 0. 
While localized heating likely is important, our results 



cannot be explained with uniform heating. Time depen- 
dent Ginzburg-Landau theory for gapped superconduc- 
tors might give an adequate phenomcnological model, 
however it is not microscopically justified in the experi- 
mental temperature range. Although our experiment is 
conceptually rather simple, a quantitative explanation is 
intriguingly complex and probably requires the incorpo- 
ration of QP diffusion. Such a model may also provide 
insight in localized nonequilibrium effects on the super- 
conducting dynamics in other cases, such as vortex mo- 
tion or flux avalanches. 
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